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Light top superpartners play a key role in stabilizing the electroweak scale in supersymmetric 
theories. For A-parity conserved supersymmetric models, traditional searches are not sensitive to 
the compressed regions. In this paper, we propose a new method targeting this region, with stop 
and neutralino mass splitting ranging from mg — m x > mt to about 20 GeV. In particular, we focus 
on the signal process in which a pair of stops are produced in association with a hard jet, and define 
a new observable Rm whose distribution has a peak in this compressed region. The position of the 
peak is closely correlated with mg. We show that for 13 TeV LHC with a luminosity of 3000 fb _ , 
this analysis can extend the reach of stop in the compressed region to mg around 800 GeV. 


Introduction With the discovery of the Higgs mm, 
particle physics reached an important milestone. How¬ 
ever, the mechanism of stabilizing the electroweak scale 
from large quantum corrections is one of the outstanding 
mysteries. In most of the models addressing this prob¬ 
lem, a key ingredient is a light top partner. As the most 
prominent example, in supersymmetry, the stop t should 
be less than about a TeV to be an effective solution to the 
fine-tuning problem [3] . Traditional searches for stops fo¬ 
cus on pair production of stops with each of them decays 
into t and the LSP, X- If m t ^ m x + m t, the top quark 
can be quite energetic. In the top pair production, a main 
background of this search, most of the top quarks are 
produced near the threshold. Based on this observation, 
various kinematical variables (e.g. mj '2 mm, Hr ,3, 
the razor variables HDHH1 and the variables invented in 
Ref. ]pG3) have been defined to distinguish stop pair pro¬ 
duction from top pair production. For recent global stud¬ 
ies of the minimal version of the supersymmetric stan¬ 
dard model, see [14J and references therein. However, in 
the compressed region where mg ~ mt + m x , the kine¬ 
matics of the top quarks from stop decay are similar to 
those in the top pair production, and such observables 
are less sensitive. In the region that m x -C mg ss mt, 
spin correlations of the top quarks can help to distin¬ 
guish the signal from background [HffiZ]- Such analysis 
has been done by the CDF, DO, ATLAS and CMS collab¬ 
orations jl8l(25j . However, with larger mg this method 
does not work well due to smaller production rate. In an¬ 
other extreme regions of the parameter space mg ~ m x , 
t decays into 4 body final states or a light quark plus 
the LSP through flavor-changing processes. In the case 
that the flavor-changing processes are important, charm 
tagging can be useful [2fil[28| . However, the jets from 
the decay are usually soft and cannot be identified. The 
leading search channel is mono-jet + MET piHTI(T| . Light 
stops can also be probed directly by comparing the ob¬ 
served tt or W + W~ pair production rate with theoretical 
calculations |4TH43] , However, it will be difficult for this 
method to be benefited from larger luminosity and higher 
energies in the future runs of LHC, since its sensitivity 


is mainly limited by systematic errors. Vector boson fu¬ 
sion tagging has also been proposed to search for stop 
in the compressed region, and it has been shown that it 
is still cannot fully close the gap in the compressed re¬ 
gion gn • II Ih e life-time of the stop is long enough, a pair 
of stops can form a bound state, the stoponium. In this 
case, searches of the stoponium can be sensitive to these 
compressed regions piol - FTTI . For recent detailed studies 
of LHC sensitivities see [JS] and references therein. In 
this region, the stops can also hadronize first and then 
decay with displaced vertices (33 ■ If the heavier stop is 
reachable, one can also study the decay of the heavier 
stop to the Higgs boson or Z 0 together with the lighter 
stop [50H53] . If the sleptons or charginos are lighter than 
the compressed stop, the decay pattern of the stop can 
be changed dramatically. See Ref. [Mj for recent study of 
these scenarios. Constraints on masses of the two stops 
can also be inferred from the measurement of the Higgs 
mass and production rate [53, 56]. Light stops also get 
constraints from low energy precision experiments such 
as the b — > sy experiment [5TH551 . 

However, it is still difficult for current searches to 
cover the compressed region with mass splitting rang¬ 
ing from mg — m x > m* to much smaller values about 
20 GeV. In this letter, we introduce a new kinemati¬ 
cal observable which targets the kinematics of this com¬ 
pressed region, and demonstrate its effectiveness. We 
note that there are other studies focusing on the similar 
parameter region [3511601 El]. In particular, the strategy 
adopted in Ref. [33] can cover parameter space around 
mg ~ mw + mb + m x , although it is less effective for 
mg ~ mt + m x . 

Kinematics around the compressed region In the com¬ 
pressed region mg > m t + m x , the t first decays into a 
pair of t and Due to the compressed nature, in the 
rest frame of t. the % and t are almost at rest. Therefore, 
in the lab frame, the transverse momenta of the t and 
X have a simple relation that pr{x) — ( m \/ m t)PT(t)- 
Therefore, in the process of the t pair production, the 
contribution to MET from the two y’s approximately 
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cancel each other, and as a result the kinematics of the 
top quarks from t pair production is very similar to those 
from the top pair production, making the search very dif¬ 
ficult. We propose to focus on events with an additional 
hard jet from initial state radiation (ISR). In this case, 
we have p T (j isr) - -(pr(ti) + Pr { i 2)) where h and t 2 
are the two stops produced in this process. Therefore, 
the ratio between MET and the pt(jisr) 


level is still around the value given by Eq. ([ 3 ]), with re¬ 
placement Am —>• |Am|. 

Similarly, we also expect to see a sharp peak in the Rm 
distribution at the compressed region toj ~ mw + m f, + 
m x , where the W boson and the LSP is approximately 
stationary in the stop rest frame. The width of the peak 
in the Rm distribution at parton level can be estimated 
using Eq. Q with nit in Eq. ^ replaced by mw- 


Rm —$t/Pt (j'isr) ~ rn x /m i , (1) 

where j/p is the total MET in this process. Hence, we ex¬ 
pect a peak-like feature in the Rm distribution. Ref. [60] 
has also noticed the similar kinematical feature. The 
spread of this peak can come from several sources. In 
the rest frame of t, the momentum acquired by x can be 
written as 


A P x = 


[( m~ — (m t + m x ) 2 )(m~ — (m t — m x ) 2 )] 1 / 2 


2 m- t 


( 2 ) 


Therefore, in the compressed region, 

< pm.Am)'/’', (3) 

where A?n = mj— m x —mt- In most part of the parameter 
space we are interested in in this study, the boosts of the 
?s are small. Therefore 


ARparton = Ap x /p T (jlSR) (4) 

can serve as a good estimate of the width of peak in the 
Rm distribution at parton level. In practice, additional 
soft radiation and detector effects will also smear the dis¬ 
tribution of Rm- Nevertheless, as we will demonstrate, 
there is still a peak in the Rm distribution around m x /mi 
in the compressed region. 

In the compressed region where < mt + m x , the 

t decays into x and b and W through a virtual t. Ne¬ 
glecting the spin correlation between the initial and final 
states the differential decay width of t with respect to the 
invariant mass of the virtual t can be written as 

dT-t _ r f\q t ) g t 2 r t (q t ) 

dq t ~ tt (g 2 - m 2 ) 2 ’ W 

where qt is the virtual mass of the top quark. T^ '( q t ) is 
the two-body decay width of t with replacement mt —> qt, 
and T t (g i ) is the decay width of t, replacing m t —> qt- 
Eq. ([5| implies that q t prefers to be as close to mt as 
possible, with maximal value gtmax = — m x . Hence, 

even the top quarks are virtual, the x decayed from each t 
still prefers to be at rest in the rest frame of t. Therefore, 
the relation shown in Eq. 0 still holds approximately. 
Since we are not far away from the region where the top 
is on-sliell, we expect the spread of the peak at parton 


SM background and basic cuts For this analysis, it is 
crucial to identify which jet is from ISR. is with a mass 
of several hundred GeV are not usually highly boosted. 
Therefore, in the compressed region, the ts from the t 
decay are also not highly boosted. Thus, the pt of the 
hardest jet in the decay chain of the t is around mt- As 
a consequence, if we require Pr(jo) mt, where jo is 
the hardest jet, we found that it is very probable that 
Jo is jisR • Hence, we will use the ratio ^t/ptUo) as an 
approximation for Rm- The requirement of a large pr(jo) 
also helps reduce the QCD background and sharpen the 
peak of the Rm distribution as shown in Eq. Q. In 
practice, we require pr(jo) > 700 GeV. 

The leptonic decay of t is always accompanied with 
neutrinos, which smears the peak structure in the Rm 
distribution. In this analysis, we focus on the hadronic 
decays and veto events with charged leptons. At par- 
ton level, 6 soft jets from top decay appears. In prac¬ 
tice, some of these soft jets may merge into a harder one. 
Therefore, we require at least three sub-leading jets with 
p T > 60 GeV. 

An important kinematical feature of the signal is that, 
in the compressed region, the pt(jisr) is approximately 
in the opposite direction as the jJy. Therefore, we require 
\ 4 >( jo ) — </>met ^7r| < 0.15. At the same time, we expect a 
significant QCD background from the mis-measurement 
of jet energy. To reduce the background due to the mis- 
measurement of the subleading jets, we require |</>met — 
4>j > 0.2 for all the jets with p T > 60 GeV. Requiring 
Pt > 60 GeV also helps to reduce the pile-up effects 
which are significant during the high luminosity runs of 
the LHC. 

To further reduce the QCD background, we require 
at least one b-jet appears in the final state. The fa- 
tagging efficiency we use is the tight b-tagging in the 
PGS detector simulator §)3], in which the b-tagging effi¬ 
ciency is about 40% for pT{b) around 100 GeV and within 
r; < 1.2. For such a tagging efficiency the mistag rate 
of light partons for CMS detector can be as small as 
0.1% [B2]. Detailed simulation shows that the main QCD 
background is from the processes with b and c quarks in 
the final states. Since it is very unlikely for b jet to be 
the leading ISR jet in the signal, we veto events with 
the leading jet passed the loose b-tagging in the PGS de¬ 
tector simulation [65] , which is around 45% for jet with 
Pt > 700 GeV. This is smaller than the current bench¬ 
mark b-tagging efficiency of the CMS and ATLAS detec- 
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FIG. 1: Contributions to background from various processes, 
after the basic cuts described in the text. 


tors, therefore our result is conservative. 

Numerical results For both the background and sig¬ 
nal, the parton level simulations are done using Mad- 
Graph5/MadEvent [63] followed by parton shower with 
PYTHIA6.4 [53]. The detector simulation is done using 
PGS4 [55] with anti-fcr jet algorithm with a distance pa¬ 
rameter of 0.5 [55]. For the background, the MLM match¬ 
ing scheme is also used to avoid double counting [67]. For 
signal we checked the results from simulations with and 
without matching that the difference is within 20%. With 
all the basic cuts discussed above, the Rm distribution 
from SM processes with the cuts previously described is 
shown in Fig. [T] The dominant contribution to the back¬ 
ground is from tt pair production with a hard ISR jet. In 
our signal region with large Rm, a significant amount of 
MET is required. Since we veto events with charged lep¬ 
tons in the final state, the dominant contribution to the 
background is from leptonic decays of top with rs, or mis- 
tagged e/fj,s. The second leading background comes from 
QCD multi-jet production with at least one of the jets 
containing a bottom or charm quark. The background 
from electroweak processes is not important due to their 
smaller rates. 

From Fig. [I] one can see that both the tt + j'isr back¬ 
ground and the QCD background exponentially decrease 
with Rm due to the lack of the source of MET. The back¬ 
ground from electroweak processes is relatively flat but 
with a suppressed rate as shown in Fig. [T] The total 
background is well fitted by a function 




FIG. 2: Rm distributions for m,- = 350 GeV on both sides of 
the m t - = m x + mt (top) and m ( - = m x + mw + rat, (bottom). 


where A and B depends on the details of the cuts, and 
in the current choice A = 47 fb and B = 5.6. 

For the signal, the Rm distribution for mi = 350 GeV 
and several different m x are shown in Fig. [2] To make 
the feature easier to visualize, we choose points very close 
to the mass thresholds with Am « 2 GeV. One can see 
that in all cases the Rm distribution is peaked at around 
m x /mi , with widths around 0.2. From Eq. |5]), the width 
generated by the phase space of the decay of t is about 
0.05. Therefore, the typical width of the peak of the Rm 
distribution induced by parton shower and detector effect 
is about 0.2. 

In order to take advantage of the peak in the Rm 
distribution, and the fact that the background decays 
exponentially with Rm we add another cut that for 
< m t + m x 


"be 

mi 


- 0.05 < 


]/>t 

Pt ( jo ) 


< 



0.15 , 


(7) 


and for mi > mt + m x 



-0.05 < 


PtUo) 


< 


mi — mt 


+0.15. (8) 


da 

dRM 


Aexp(—BR M ) , 


As the m x —> 0, the background in the cut window de¬ 
fined in Eq. ([T]) grows exponentially. Therefore, for the 



























































































































4 


800 


600 


> 

a) 

2 400 

\ 

£ 

200 


0 

200 400 600 800 1000 200 400 600 800 1000 

m~ (GeV) m~ (GeV) 

FIG. 3: The projected 5 a discovery reach (left) and 95% C.L. exclusion limit (right) of stop in the compressed region. 




region mi > mt + m x , we choose a different window as 
shown in Eq. (|8j) which is independent of m x , so that 
more parameter space in the bulk region can be covered 
effectively. We define the 5cr and 2cr expected limit by 
S/VB = 5 and S/^/B + S = 2, where S and B are the 
signal and background after all the cuts. The 5 a reach 
at CMS with the center-of-mass energy of 13 TeV and 
luminosities of 300 and 3000 fb^ 1 are shown in left panel 
of Fig. [3] together with the expected 5er sensitivity of the 
direct stop pair production at 14 TeV LHC with CMS 
detector [BSC The right panel shows the 2er expected 
95% C.L. exclusion limit together with the current com¬ 
bined limit from direct stop pair production I32IE8EI]. 
The current limits and prospective reach of ATLAS are 
similar to CMS [25} I72H75] , which are not shown here. 
One can see that t in the compressed region with a mass 
around 600 GeV can be discovered in the LHC with 3000 
fb . It can also exclude stop with mass up to about 800 
GeV. Notice that there are ’’spike-like” features around 
the thresholds, mj — m x ~ mt and toj — m x ~ mw + n%b, 
due to the fact that the peak in Rm is sharper around 
these thresholds. 

Conclusion and discussion We point out a useful 
kinematical feature in the production of stop associate 
with an ISR jet, which can enhance the sensitivity in 
the compressed regions, with mass splitting ranging from 
mi — m x ~ mt to about 20 GeV. We show that in this 
region the observable Rm defined in Eq. ([!]) has a peak 
around Using this kinematical feature, we es¬ 

timated that this gap can be covered up to around 800 
GeV with 13 TeV LHC at a luminosity of 3000 fb -1 . 


Although we have focused on the stop searches, the 
same technique is obviously applicable to the search of 
other top partner signals with similar final states. 

In the discussion, we neglect the flavor changing decay 
mode t —>• yc, since it is model dependent and strongly 
constrained by flavor physics. It has been shown that in 
the minimal flavor violation scenario the branching ratio 
of this process is often subdominant to the four body 
decay of the stop in the region that < m x + mw + 
m;, m. For detailed next-to-leading-order studies of the 
stop decay pattern in this region, see HU EH III]- This 
method is also applicable to other decay chains of stop 
in the compressed region EH}- 

The main background of this analysis is from top pair 
production associated with ISR jets, with at least one of 
the top decays leptonically and the charged lepton fails 
the lepton veto. A majority of these events have the 
charged lepton close to hadronic activities. One may be 
able to distinguish these events further with alternative 
lepton isolation criteria. The top quarks generated from 
the stop decay in the signal are in general with a smaller 
boost than the top quarks from the background since 
mi > mt , which is assumed in this analysis. We may be 
able to use this property to further distinguish the top 
quarks from the signal and from the background. We 
leave these analysis to future work. 
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